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Abstract 

A hybrid formalism is proposed for interacting classical and quantum sytems. This 
formalism is mathematically consistent and reduces to standard classical and quantum 
mechanics in the case of no interaction. However, in the presence of interaction, the 
correspondence principle is violated. 
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Quantum mechanics gives exceedingly accurate predictions for atomic and nuclear sys- 
tems. Classical mechanics is just as successful for planetary motion. Can there be a 
theory encompassing both classical and quantum mechanics, so that atoms and planets 
can be treated in a unified way? Indeed, it was speculated long ago that the gravitational 
field is classical, even though its material sources are quantized pj. Niels Bohr always 
insisted that measuring instruments must be described in classical terms f2|. However, 
Bohr did not provide a dynamical description of their interaction with quantum systems. 
A dual (classical and quantum) description of measuring instruments is possible J5|, but 
this is not a truly hybrid formalism as proposed below, because in that dual formalism 
the interaction is always between subsystems of the same kind. 

If we try to write equations of motion that combine classical canonical variables and 
quantum operators, we find that the classical variables, whose values ought to be ordinary 
numbers, pick up operator components, which is manifestly unacceptable (we shall later 
show how to overcome this difficulty). There have been numerous claims and counter- 
claims as to the possibility of existence a consistent classical-quantum formalism. Some 
recent articles are listed in the references [Q . There is no real contradiction between them, 
because their conflicting conclusions result from different demands for consistency. 

In this article, we propose a hybrid formalism for interacting classical and quantum 
sytems. This formalism is mathematically consistent and reduces to standard classical 
and quantum mechanics if there is no interaction between the systems. However, in 
the presence of interaction, the correspondence principle is violated, so that this hybrid 
formalism appears physically abnormal. 

Our work follows a remark by Jauch || that if we assume that not all quantum dy- 
namical variables are actually observable, and if we set rules for distinguishing measurable 
from nonmeasurable operators, it is then possible to define a classical system as a special 
type of quantum system for which all measurable operators commute. All the observable 
dynamical variables of such a classical system can be simultaneously diagonalized and 
therefore can, in principle, have sharp values. If we are confined to this set of observables 
and we use a basis where they all are diagonal, then only diagonal elements of density ma- 
trices can be reconstructed from measurement statistics and have an operational meaning. 
Quantum states (i.e., density matrices) with the same diagonal elements are indistinguish- 
able. We may therefore consider all states as pure by defining their off-diagonal elements 
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as p mn = y/PmmPnn, so that p always has rank 1. It is then possible to define a classical 
state vector ip by p = tpipK 

Pure states defined in this way do not imply that dynamical variables have sharp 
values. The situation is just as in classical statistical mechanics: the expectation value of 
an observable A is given by 

(A) =tr (pA), (1) 

and we may then have (A 2 ) > (A) 2 , as usual. Note that such a "classical" system also has 
noncommuting operators, but the latter should be considered as abstract mathematical 
expressions which are not experimentally observable. With this formal definition, the 
notion "classical" acquires a meaning with respect to a specified set of dynamical variables. 
This does not preclude that, by using more elaborate observation means, every physical 
system might display quantum features. 

To completely mimic classical mechanics, we need a quantum algorithm that reproduces 
exactly the equations of motion of a classical canonical system. For this, we shall follow 
Koopman's formalism || [jj. We consider for simplicity a single degree of freedom and 
denote the canonical variables as x and y (rather than x and p, as usual, because we wish 
to reserve the symbol p for the momentum of a quantum system, to be introduced later). 
Let us write the Liouville equation as 

idf/dt = Lf, (2) 

where L is the Liouville operator, or Liouvillian, 

L = (-^) - (H) {-%) ■ (3) 

The Liouville density / is never negative. Since its quantum analog, namely Wigner's 
distribution JF|, [8| (which may be negative), is quadratic in the quantum wave function, 
it is convenient to introduce likewise a "classical wave function" ip = v7, which in this 
case satisfies the same equation of motion as /, 

idi/j/dt = Lip. (4) 

We shall now consider ip as the fundamental entity (but only / = \ip\ 2 has a di- 
rect physical meaning). It can be proved that, under reasonable assumptions about the 
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Hamiltonian, the Liouvillian is an essentially self-adjoint operator and generates a unitary 
evolution |§: 

(0,0):= J ip(x, y, t)* 0(x, y, t) dxdy = const. (5) 

The introduction of a Hilbert space with unitary dynamics enables one to use familiar 
methods of quantum mechanics for the analysis of ergodic problems [U and of classical 
chaos [HJ. It is possible to further mimic quantum theory by introducing commuting 



operators x and y, defined by 

xif) = xi/;(x, y, t) and yip = yip(x,y,t). (6) 

Note that the momentum y is not the shift operator (the latter is p x = —id/dx). Likewise 
the boost operator is p y = —id/dy. These two operators are not observable. We shall 
henceforth omit the hats over the classical operators, as there is no danger of confusion. 

The analogy with quantum mechanics can be pushed further. What we have above is 
a "Schrodinger picture" (operators are constant, wave functions evolve in time as if)(t) = 
U(t)ip(0), where U(t) = e~ tLt if the Hamiltonian is time-independent). We can also define 
a "Heisenberg picture" where wave functions are fixed and operators evolve: 

X H (t) = U^XU. (7) 

The Heisenberg equation of motion, 

i dX H /dt = [X H , L H ] = U ] [X, L] U, (8) 



together with the Liouvillian (S), readily give Hamilton's equations 
dx dH dy OH 



(9) 



dt dy ' dt dx 

There is however an important difference: the time translation operator L is not the 
energy, and its spectrum may extend to — oo. For example, if we have a harmonic oscillator 
with H = (x 2 + y 2 )/2, the Liouvillian is 

L = yp x - xp y , (10) 

whose eigenvalues are all the integers, negative as well as positive. There is nothing wrong 
in that, since L involves the unobservable shift operators p x = —id/dx and p y = —id/dy, 



and therefore L itself is not observable. Note that the solution of Hamilton's equations 
does not introduce non-observable components into the observable variables, because the 
Hamiltonian (contrary to the Liouvillian) involves only the observables x and y. 

It is easy to introduce into the above dynamical formalism a quantum system, with 
conjugate dynamical variables q and p, as long as the two systems do not interact. In the 
Schrodinger picture, we have a wave function ip(x,y,q) whose evolution is given by 

dip/dt = Kijj, (11) 

where K is "Koopmanian" 

K — L(x, y,p x ,P y ) + H(q,p). (12) 

In that wave function, the coordinates x, y, and q may be entangled. In any case, we can 
obtain the reduced density matrices of the "classical" and quantum systems by means of 
partial traces on the other system, as usual. 

The non-trivial problem is to introduce an interaction between the two systems. An 
apparent difficulty occurs in the Heisenberg picture because the equations of motion will 
mix variables of all kinds: classical observables, non-observables, and quantum operators. 
This still is acceptable, because we do not predict actual values for these variables, but 
only expectation values, (A) = tr(pA), and no contradiction may occur. Anyway, we 
can always use the Schrodinger picture, where we know how to handle entangled wave 
functions. 

The true difficulty, as we shall now see in a simple example, is that the correspondence 
principle fails. Consider two harmonic oscillators, with a bilinear coupling kqx, where k 
is a constant. If we treat both of them classically, with a Hamiltonian 

H = l(q 2 +p 2 + x 2 + y 2 ) + kqx, (13) 

we obtain equations of motion 

q=p, p=-q-kx, (14) 

x — y, y — —x — kq, (15) 

and there are two characteristic frequencies, to = y/l ± k, corresponding to the normal 
modes (q±x). Exactly the same equations of motion appear in the Heisenberg picture for 



quantum mechanics. Likewise, in the Schrodinger picture, whenever there is a quadratic 
potential, the differential equation for the Wigner distribution is identical to the Liouville 
equation in classical mechanics J/J. In view of this formal agreement of classical and 
quantum mechanics, it is natural to demand that any hybridization of the system in 
Eq. (pT3|) shall also give the same equations of motion. This is the definite benchmark we 
propose for an acceptable classical- quantum hybrid formalism. It is not obvious that this 
criterion can be achieved, because the hybrid system has no normal modes (g ± x). 
Let us try to obtain the above equations of motion from a Koopmanian 

K = \ (q 2 +p 2 ) + (yp x -xp y ) + K i , (16) 

with a suitably chosen interaction term Ki. We cannot have both \p, Ki] = —kx and 
[y, Kj\ = —kq together with [y,p] = 0, because these equations are incompatible with 
Jacobi's identity 

[y, \p,Ki]] + \p, [K it y]] + [K h [y,p]} = 0. (17) 

The best result we were able to obtain was by "koopmanizing" the interaction term kqx 
in the Hamiltonian (13). We have, from Eq. (^), 



Ki = ikqd/dy = —kqp y . (18) 

The resulting equations of motion are Eq. (p~5|) , unchanged, and 

q = p, p=- q - kp y , (19) 

Px = Py, P y = ~Px- (20) 

The last equation is necessary, because an unobservable variable p y appears in Eq. ([19]) • 
The solution of (20) for p y is a superposition of sint and cost terms. When substituted 
into Eq. (|T9|), these terms behave as a driving force with resonant frequency, so that q and 
p include terms behaving as tsmt and tcost. The amplitude of the quantum oscillator 
increases linearly with time, and energy is not conserved. It need not be: what is conserved 
is the Koopmanian, which includes a term (yp x — xp y ) whose spectrum extends to — oo. 

In conclusion, there is no mathematical inconsistency in the hybrid formalism that 
we proposed. However, it violates the correspondence principle which we would expect 
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to hold exactly for a pair of oscillators with bilinear coupling. Therefore such a theory 
appears quite abnormal from the point of view of physics. 

We thank Lajos Diosi for stimulating discussions. DRT was supported by a grant from 
the Technion Graduate School. Work by AP was supported by the Gerard Swope Fund 
and the Fund for Encouragement of Research. 
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